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Abstract 

In this article we prove the existence of solutions to the singular coagulation equation with multifrag- 
mentation. We use weighted //-spaces to deal with the singularities. The Smoluchowski kernel is covered 
by our proof. The weak L 1 compactness methods are applied to suitably chosen approximating equations 
as a base of our proof. A more restrictive uniqueness result is also given. 



1 Introduction 

The coagulation process describes the kinetics of particle growth where particles can coagulate to form larger 
particles via binary interaction. On the oder side, the fragmentation process describes how particles break into 



two or more fragments. Examples of these processes can be found e.g. in astrophysics, polymer science 15 
and cloud physics (lOj Chapter 15]. 

The dynamic of the coagulation-fragmentation process is described by the integro-differential equation 

X oo 

du(x,t) If f 

— 7^— = 2 / K {x-y,y)u{x-y,t)u{y,t)dy- \ K(x i y)u(x i t)u(y,t)dy 
o o 

b(x, y)S(y)u(y, t) dy - S(x)u(x, t), (1) 

X 

with initial condition 

u(x,0) = u Q {x) > a.e. (2) 

where the non-negative variables x and t represent the size of the particles and time respectively. The values 
u(x, t) denote the number density of particles with size x at time t. The rate at which particles of size x coalesce 
with particles of size y is represented by the coagulation kernel K(x,y). The rate at which particules of size 
x are selected to break is determined by the selection function S(x). The breakage function b(x,y) gives the 
number of particules of size x produced when a particule of size y breaks up. 

The equation ([lj is named the continuous coagulation equation with multiple-fragmentation where the 
multiple-fragmentation kernel T defines the selection function S and the breakage function b by 

X 

S(x) = -T(x,y)dy, b(x,y) = T(y,x)/ S(y), 

J y 



or vice versa. 

The breakage function is assumed here to have the following properties 

y 

xb(x,y)dx = y for all y > 0, (3) 

o 



* Corresponding author. Fax: +49 391 6718073, 
E-mail address: cuetacam@st.ovgu.de (C. Cueto Camejo). 



1 



which is the conservation of mass and 
y 

J b(x, y)dx = N <00 for all y > 0, b(x, y) = for x > y, (4) 
o 

where the parameter N represents the number of particles produced in fragmentation events. In this paper N 
is assumed to be finite and independent of y. Equation (|3| allows the system to conserve the total mass during 
the fragmentation events. It states that the total mass of the fragments is equal to the mass y of the particle 
that breaks. 

The existence and uniqueness of solutions to the coagulation-fragmentation equation has already been the 
content of several papers results. Nevertheless, the case of multifragmentation, that is, when the particules can 
break into two or more parts, has not yet been treated very often. Most results concern binary breakage. The first 
result in regard to the existence and uniqueness of solutions to the coagulation equation with multifragmentation 
was given by Melzak (8j, where the coagulation and multifragmentation kernels were assumed to be bounded. 
For more recently result see e.g. [7], |6j and [5j. In [7], Laurengot considered unbounded coagulation and 
fragmentation kernels, for coagulation kernels on the form K(x,y) = r{x)r{y) with no growth restriction on 
r and a small growth asumptions on T depending on r. Giri et al. |5j studied the coagulation kernels of 
the form K(x,y) — (f>(x)(f)(y) for some sublinear function <fi under the growth restriction cj>(x) < (1 + x) 11 for 
< fx < 1, and the selection function S(x) is also considered under the same growth assumption. In [6j, 
Giri et al. proved the existence of solutions to the coagulation equation with multifragmentation for a more 
general fragmentation kernel, in order to cover the fragmentation kernel T(y,x) = (a + 2)x a y 1 ~^ a+1 ' > getting 
a result for a > — 1 and 7 G]0, a + 2[. The existence proofs in [6j and |5j are based on the well known basic 
method by Stewart |12| , where the solution is obtained through the convergence of the solutions to a sequence 
of truncated problems. In [6] a nd [5j the existences and uniqueness of solutions to the truncated problems 
is claimed to be solved by |12[ Theorem 3.1], but the truncated problems in [6] and |5j do not satisfy the 
hypotheses of |12[ Theorem 3.1], leaving in that way a gap in the proof of the original problems. Stewart 
required in (l2{ Theorem 3.1] the binary fragmentation kernels F(x,y) to be symmetric, i.e. F(x,y) = F{y,x) 
while the multifragmentation kernels T(y,x) in (HI and (HI are not necessarily symmetric. On the oder hand, 
in [6] and [5] the study is concerned on the multifragmentation kernels which is a more general case than the 
binary fragmentation studied by Stewart [12]. In (6j the uniqueness of the solutions was not studied. 

In the present article, our first aim is to prove the existence and uniqueness of solutions to the coagulation 
equation with multifragmentation with singular coagulation kernels 

K(x,y) <k{l + x) x {l + y) x (xy)- a with A - a G [0, 1[, a G [0, 1/2], (5) 

giving in this way an existence and uniqueness result for the case of the important Smoluchowki coagulation 
kernel 

K(x, y) = (x 1 / 3 + y 1 / 3 )^- 1 / 3 + jT 1 / 3 ) 
for Brownian motion, see Smoluchowski |11] , The equi-partition of kinetic energy (EKE) kernel 



K{x,y) = {x^ + y^f^\ 



•I 1 

x y 

is also covered by our analysis. We want to point out that there is no previous result concerning the coagulation- 
fragmentation equation with singular coagulation kernels. Our existence result is based on the proof of Stewart 
[12] . We extend the methods we developed in [2] for singular kernels in the pure coagulation problem. 
For our existence and uniqueness result we consider the class of fragmentation kernels 

v 

r(y, x) < y e b(x, y) with J b(x, y)x- 2a < Cy' 2 * for 9 G [0, l[,cr G [0, 1/2], and a constant C. (6) 



The class of frangmentation kernels ^ includes the kernel T(y, x) — (a + 2)x a y 1 ~(- a+1 '> for a > — 1 and 7 G [0, 1[. 
This kernel was studied by Giri et al. |6j, where they proved the existence of weak solutions to the coagulation 
equation with multifragmentation, but with nonsingular coagulation kernels. 

Our second aim is to fill the gap left in (6j. For a = kernels ([5j and (|6j become 

K(x,y) < k(l + x) x (l + y) A with Ag[0,1[ (7) 

and 

T(y,x)<y e b(x,y) with 0G[O,1[. (8) 
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The coagulation kernels ([7]) are the ones considered in |6] and taking b(x, y) = x~ T in Q with r G [0, 1[ we also 
have the fragmentation kernel from [6] . Therefore our proof of existence and uniqueness of the solutions to the 
truncated problem in |6] will be covered by our result for a — 0. 

In order to study the existence of solutions of ([lJ-Q. we define for some given a g [0, 1/2] the space Y to 
be the following Banach space with norm || ■ ||y 



oo 

Y = {u€ £ 1 (]0, oof) : < oo} where ||u|| y = J (x + x~ 2tr )\u{x, t)\dx. 



That Y is a Banach space is easily seen. We also write 



(i 



\\u\\ x = j xu(x,t)dx and UuU^-an — J x a u(x,t)dx, 



and set 

Y + = {ueY:u>0 a.e.}. 

Now we define a weak solution to problem ([l])-(|2| in the s ame way as Stewart fl2] : 

Definition 1.1 Let T G]0, oo]. A solution u(x,t) of ^jty-^ty is a function u : [0, T[ — > Y + such that for a.e. 
x G [0, oof and t G [0, T[ the following properties hold 



i) u(x, t) > for all t G [0, oof, 

ii) u(x, •) is continuous on [0, T[, 

iii) for all t G [0, T[ the following integral is bounded 



K (x, y)u{y, r) dy dr < oo and J J b{x, y) S '(y)u(y, r) dy dr < oo 

Ox 

iv) for all t G [0,T[, u satisfies the following weak formulation of |7p 



u(x, t) = u(x, 0) + 



o 



K(x - y,y)u(x - y,T)u(y,T) dy - / K{x,y)u{x,T)u[y,r) dy 

( 



WK*,»)S(»MM)*-*MM) - 



In the next sections we make use of the following hypotheses 
Hypotheses 1.2 

HI) K{x,y) is a continuous non-negative function on ]0, oo[x]0, oof, 

H2) K{x,y) is a symmetric function, i.e. K(x,y) = K(y,x) for all x,y G]0, oof, 

H3) K(x,y) < K{l + x) x (l + y) x (xy)- a for X - a G [0, 1[,ctG [0,1/2], and constant k, 

H4) S(x) <y e for 0G [0,1 [, 



v 

2a 



H5) b(x,y) is such that J b(x,y)x 2<T dx < Cy 

o 

In the rest of the paper we consider k = 1 for the simplicity. 

We study the uniqueness of the solutions to (JlJ— ([2J under the following further hypotheses 

H3') K(x, y) < K 1 (x- a + x x - a )(y- a + y x ~ a ) such that a, A - a G [0, 1/2] and ki > 0, 
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H4') S(x) =y e for 0<X-a. 

The restriction A — a € [0, 1/2] in H3') limits our uniqueness result to a subset of the kernels of the class 
defined in H3), namely to the ones for which A — a E [0, 1/2] holds. But the class of kernels defined in H3') 
is also wider than the defined in H3) for A — a <E [0, 1/2]. In this way we are also giving uniqueness result for 
kernels which are not included in the class defined in H3). On the oder hand, the restriction 9 < A — a in H4') 
limits our uniqueness result to a more restricted class of fragmentation kernels. 

We introduce now some easily derived inequalities that will be used throughout the paper. The proof of 
these inequalities can be found in Giri |4j. For any x, y > 

2 P~ 1 ( X P +yP) < ( X + y)P < X P + yP if < p < 1, (9) 

2 p ~ 1 (x p + y p )> {x + y) p >x p + y p if p > 1, (10) 
2 p -\x p + y p ) > (x + y) p if p < 0. (11) 

The paper is organized as follows. In Section 2 we define the sequence of truncated problems and prove in 
Theorem |2.4| the existence and uniqueness of solutions to them. We extract a weakly convergent subsequence 
in L 1 from a sequence of unique solutions for truncated equations to ([lJ-Q. In Section 3 we show that the 
solution of (JT|) is actually the limit function obtained from the weakly convergent subsequence of solutions of the 



truncated problem. In Section 4 we prove the uniqueness, based on the method of Stewart 13 , of the solutions 
to ([lJ-([2]) for a modification of the classes §5§ and ^ of coagulation and fragmentation kernels respectively. We 
obtain uniqueness for some kernels which are not covered by the existence result. 



2 The Truncated Problem 



We prove the existence of a solution to the problem ([lJ-Q by taking the limit of the sequence of solutions of 
the equations given by replacing the kernel K(x,y) and the selection function S(x) by their respective 'cut-off' 
kernel K n (x,y) and S n (x) for any given n G N 



K n (x,y) = 



K(x,y) if x + y<n and x,y>a/n 
otherwise. 



S n {x) 



S(x) if x < n 
otherwise. 



(12) 



Note that if we take a = our truncated problem will be defined as in Giri et al. j6J. For the defined kernels 
the resulting equations are written as 



du n (x,t) 
dt 



i J K n (x - y 7 y)u n (x - y,t)u n (y,t) dy - J K n (x,y)u n (x,t)u n (y,t) dy 
o o 

n 

+ f S n {y)b(x,y)u n (y,t)dy - S n (x)u n (x,t), 



with the truncated initial data 

<{x) - 



Uq(x) if < x < n 
otherwise, 



(13) 



(14) 



where u n denotes the solution of the problem ( 13 1-( 14 ) for x £ [0, n]. Next, we rewrite our truncated problem 
(13l-(14l in an equivalent form. We prove some lemmas, which are used to show the existence and uniqueness 
of the solution to the truncated problem. 

2.1 Existence and uniqueness of solutions of the truncated problem 

Let us define the operator P as 



P{x,t,u n ) 



K n (x, y)u n (y, r) dy + S n (x) 



dr. 
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which allows us to rewrite the truncated problem ( 13 )-( 14 1 in the equivalent form 

K n (x - y)u n (x - y, t)u n (y, t) dy 



~ [u n (x,t)exp(P(x,t,u n ))] = ~exp(P(x,t,u n )) 



S n {y)b(x,y)u n (y,t)dy 



with 

<{x) - 

Now, we define the operator G as 



uq(x) if < x < n 
otherwise, 



(15) 



(16) 



G(c)(x, t) = i / exp ( - [P(x, t, c) - P(x, r, c)] ) 



K n {x - y, y)c(x - y, r)c(y, r) dy 



s n (y)H x iy) c (y, T ) d v 



dr + Uq(x) exp ( — P(x, t 7 c)) , 



(17) 



for c € C ([0, T]; L 1 (]0, n[)j . By using (151 and ( IT I , it can be easily checked that a solution u n to (13l-(14l 
satisfies 

u n {x,t) = G(u n )(x,t). (18) 

The problems ( 13 l -(|l4| ) and ( 16 )-( 18 1 are equivalent. As a consequence, we prove the existence of solutions of 
the problem ( 16 1-( 18 ) . With this aim, we use the contraction mapping principle in some interval [0,T]. But 
first, we introduce some necessary definitions. Let us set 

M = max {sup {X„ (a;, y) :i,i/G [0, n}} , sup {S n (x)b(x, y) : x,y G [0,n]}} , (19) 
L = {n a <j-°MT + X)\\v%\\y, 

and choose ti,t 2 > such that 

exp (2n 2rT MLt 1 ) (l + 2 1 - 2a n 6a a- 2 ' 7 MLh + 2Cn e t 1 ) < 2 (20) 
t 2 exp (2n 2a MLt 2 ) [2 1 - 2a n 6 ' 7 cr" 2<T ML(Lt 2 + 1) + Cn e (2Lt 2 + 1) + n 2a ML] < 1. (21) 

We set 

t = min(ti, t 2 , T). 

l 

For those c G C([0, T]; i 1 Q0, n[)) for which / x~ 2a \c{x, t)\dx is finite for all t e [0,t }, we define the norm || • \\n 

o 

by 

n 

\\c\\d — sup / x~ 2cr \c(x, t)\ dx. 
te[o,t ] J 



Now we set 

D = {ce C([0,t ];£ 1 aO,nD) : ||c|| D < 2L] . 
Then, by defnition of P and the non-negativity of for c G D we have 

/ t n—x 

exp (- [P(x,t,c) - P(x,t,c)]) = exp I - / / K n (x, y)c(y, s) dy + S(x) 



o 

t n—x 



ds 



< exp 



K n (x,y)c(y,s)dyds 



T 



< exp (n 2a M\\c\\ D t) < exp (2n 2a MLt). 
For the proof of Theorem |2.4| some lemmas are necessary which we present now. 



(22) 
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Lemma 2.1 The functional G maps the set D into itself. 



Proof. Choose c such that ||c||d < 2L. For t £ [0, fo]> using (17 1, (22 1, the definition of M, and Fubini's 
Theorem, we get 



\G(c)(x,t)\x~ 2a dx 

exp(-[P(x,t,c)-P(x,T,c)]) 

n 

+ J S n (y)b(x,y)c[y,r) dy 
K n (x,y)\c{y,s)\dy ds 
K n (x-y,y)\c(x-y,T)\\c(y,T)\dy+ I S n {y)b{x,y)\c(y,T)\dy 



Lo 



K n {x - y, y)c{x - y, r)c(y, r) dy 
dr + Uq(x) exp ( - P(x, t, c)) 



x~ 2rT dx 



< 



exp 



x 2a dx dr 



+ / K(x)|exp / / K n (x,y)\c(y,s)\dyds x 2a dx 



vO 
t f n x 



<^exp(2n 2a MLt) J J J K n (x - y,y)\c(x - y,T)\\c(y,T)\x- 2a dy dx 



o Lo o 



S n (y)b(x,y)\c(y,T)\x 2a dydx 



x 



dr + \\uq\\y exp (2n 2cr M Lt) . 



Changing the order of integration, then a change of variable x — y — z, and again re-changing the order of 
integration while replacing z by x gives 



\G(c)(x,t)\x~ 2lT dx 



< exp (2n 2a MLt) 



IWIIy + 2 



n n—x 



.0 



K n (x,y)\c(x,T)\\c(y,T)\(x + y) a dydx 

n y 

+ j l Sn ^ h ( X -i y ^ C ( y,T ^ X ~ 2a dxdV 




dr 



(23) 



By using the definition of K n and S n we have 
\G{c){x,t)\x~ 2a dx 

n—a/n n - 



< exp (2n 2a MLt 



\<\w + 



a/n (7 / Ti 



K(x,y)\c(x,T)\\c(y, t)\(x + y) 2a dydxdr 

o 

n y 

+ Jfs(u)hi.r.!,)\r(!I.T)\.r 2 "d.r<h, 





dr 



(24) 



(. 



Now, we multiply and divide by (xy) 2a inside the first integral term on right hand side of (24 1. Taking in 



account the definition of M (19) and H5) in Hypotheses 1.2 we get 
\G(c)(x,t)\x- 2a dx 



< exp (2n 2a MLt) 



K\\y + 



2 -(l+2 CT)71 6 V -2a M 



n — a/n n—x 



\c(x,T)\\c(y,T)\(xy) 2a dydxdr 



a/ 7i a/n 



+Cn 9 J \c(y,T)\y- 2 °dy 
o 



d.T 



Now, using that c <E D we obtain 

\G(c)(x,t)\x~ 2a dx < exp(2n 2a MLt) (\\u%\\ Y + 2- ( - 1+2<T) n 6a cj- 2a Mt\\c\\ 2 D + Cn 9 t\\c\\ D 
< exp {2n 2a MLt) (l + 2 1 ~ 2<T n &a a^ 2(y M Lt + 2Cn e t) L. 



Then, by using (20 1 we find 



(25) 



\G{c){x,t)\x- 2a dx < 2L. 

Hence, by definition of || • \\n we have ||G||d < 2L and this completes the proof of the lemma. 

Lemma 2.2 Consider ci, c 2 £ C([0, to]; i 1 ]0, n[) and let B = max {\\ci\\d,\\c2\\d}, i.e. B < 2L as well as 

H(x,T,t) = exp ( - [P{x,t,c{) - P(x,t,c 1 )] ) - exp ( - [P(x,t,c 2 ) - P(x,T, c 2 )] ). 
Then, for < r < t < to and < x < n we have 



□ 



\H(x,T,t)\ < (t-r)n 2<T Mexp((t-T)n 2<T SM)||ci - c 2 || 



(26) 



Proof: For any fixed x £ [0, n], t,r £ [0, to] we may assume c\ to satisfy 

P(x, i, ci) - P(x, t, ci) > P(x, t, c 2 ) - P(x, r, c 2 ). 

Then 

|H(a:,T,t)| = -H(x,r,t) 

= exp ( - [P(x, t, c 2 ) - P(x, r, c 2 )] ) 

■ [1 - exp ( - [P(x, t, ci) - P(ar, r, ci) - [P(x, i, c 2 ) - P(x, r, c 2 )]])] . 



(27) 



Since 1 — exp ( — x) < x for x > 0, (27 1, together with the definitions of B and M, and the non-negativity of 
S n {x) leads to 

|i/(a;,T,t)| < exp ( - [P(x, t, c 2 ) - P(ar, r, c 2 )]) [P(ar, i, ci) — P(x, r, ci) — t, c 2 ) — r, c 2 )]] 
= exp (- [P(x,t,c 2 ) - P(x,t,c 2 )]) / / K n (x,y)[a(y,s) - c 2 (y : s)]dyds 



< 



exp 



K n (x,y)c 2 (y,s)dyds (t - r)n 2cr M || ci -c x || 



D 



r 
,2ct 



<(t- r)n 2a M exp ((t - t)h 2 °5M) || Cl - c 2 | 



□ 
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Lemma 2.3 For c l5 c 2 £ D and to as above there exists 76 [0, 1[ such that 

||G(c 1 )-G(c 2 )|| D < 7 ||c 1 -c 2 |b, 
i.e. the operator G is a contraction. 

Proof. Choose c±, c 2 £ D. Using the defintion of G we find that 

G( Cl )-G(c 2 ) = u%(x) [exp(-P(x,t,c 1 ))-exp(-P(x,t,c 2 ))] 

t X 

+ 2 / ex P ( ~ [P(x,t,Ci) - P(x, r, ci)] ) K n (x - y,y)c 1 (x - y,T)d(y,T)dydT 



exp ( - [P(x,t,c 2 ) - P(aj ! r,c 2 )] ) y #„(a; - y,y)c 2 (x ~ y,T)c 2 (y,T) dy dr 



n 

exp(- [P(x,t,ci) - P(x,t,Ci)}) / S n (y)b(x,y)c 1 (y,T)dydT 



cxp(- [P(x,t, c 2 ) -P(x, r, c 2 )]) / S n (y)b(x,y)c 2 (y,T)dydT 



By addition and subtraction of the terms 



exp ( - [P(x,i,c 2 ) - P(x,t,c 2 )] ) / K n (x - y,y)c\{x - y,T)ci(y,r) dy dr 



and 



exp (- [P(x,t,c 2 ) - P{x,t,c 2 )]) / S n (y)b(x,y)c 1 (y,T)dydT, 



together with the defintion of H in Lemma 2.2 we rewrite (28 1 as 

t 



G(ci) - G(c 2 ) - t#(a:)H (s, t, 0) + / JTfo t, r) 



#n(z ~ 2/> J/)ci(a; - t/, r)ci(j/, r) dy 



dr- exp ( - [P(x, t, c 2 ) - P(x, r, c 2 )] ) 



S„(y)b(x,y)c 1 (y,r) dy 



K„{x - y, y)c 1 (y, r) [c 2 (a; - j/,r) - c x (x - y, r)] 



+ 2 / K n( x -y,y)c2(x-y,T)[c 2 (y,T)-c 1 {y 1 T)]dy 



Sn(y)H x ,y) [°2{V,T) - c 1 (y,r)}dy 



dr. 
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Now, making use of the definition of || • \\ x -2a it follows that 

||G( Cl )-G( C2 )|| x _ 2( , 



Tl Z X 

J Uo(x)H(x,0,t) + J H(x,r,t) ^ J K n (x~y,y)c 1 (x-y,T)c 1 (y,r)dy 

L 

n It 

+ J S n (y)b(x,y)c 1 (y,r)dy dr - J exp (- [P(x,t,c 2 ) - P(x,t,c 2 )]) 

x Jo 

x 

* J K n (x - y,y)ci{y,T)[c2{x- y,r) - c^x - y,r)]dy 
. o 

x 



+ J S„(y)b(x,y)[c 2 (y,T) - ci(y,r)]dy 

X 

Applying the triangule inequality we obtain 

||G( Cl )-G(c 2 )|| x _ 2 . 



dr 



x~ 2a dx. 



< 



x 

n x 



i u x 

j \H(x,r,t)\ \jj K n (x - y,y)\ Cl (x - y,T)\\ci(y,T)\x~ 2a dydx 
o Loo 

n n It 

J J '-S , n (j/)6(a;,j/)|ci(j/,T)|a; _2<7 dyda; dr - J exp ( - [P(x,t, c 2 ) - P{x,t, c 2 )] ) 

o 

if„(x - y,y)\ci{y,T)\\c 2 {x -y,r)- a(x - y , r)\x^ 2a dy dx 

X 

\J J K n(x - y,y)\c 2 (x - y,T)\\c 2 (y,T) - ci(y,T)\x^ 2rT dydx 



n n n 

J J S n (y)b(x,y)\c 2 {y,T)- Cl {y 7 T)\x- 2 < T dydx dr + [ u%{x)\H(x,0,t)\x- 2a dx. 



o o 

n x 



x 
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Now, changing the order of integration, then a change of variable x — y = z, then re-changing the order of 
integration while replacing z by x we have 

\\G{c x )- G(c 2 )\\ x . 2 „ 



n n—x 



<J\H(x,T,t)\ X -J y" /f„(x,y)|c 1 (a:,T)||c 1 (y,T)|( a : + 2/ )- 2CT dy^ 







+ / / Sn{y)b{x,y)\ci(y,T)\x a dydx 



x 



z 

dr - J exp(- [P(x, t, c 2 ) - P(x, t, c 2 )] ) 



n n—x 





n n—x 





n n 



K n {x,y)\ci{y,T)\\o2{x,T) - c\ (x, r) | (x + y) 2a dydx 
K n (x,y)\c 2 (x,r)\\c 2 (y,T) - ci(y,r)|(x + y)~ 2a dy dx 

dr+ I u%(x)\H{x 7 0,t)\x~ 2 ° dx. 



+ j j S n (y)b(x,y)\c 2 {y,T) - a(y,T)\x a dydx 

x 

By using Lemma [2. 2| it gives 



\\G( Cl )-G(c 2 )\\ x - 2a 

< (i-T)n 2<T Afexp(2(t-r)n 2cr ML)||ci - c 2 \\ D 
t 



n n—x 



L 



K n (x,y)\ci(x,T)\\ci(y,T)\(x + y) a dy+ / / S n (y)b(x,y)\ci(y,T)\x a dxdy 



dr 



+ exp(2tn 2 °ML) J J K n (x,y)\c 2 (x,T)\\ Cl {y,T) - c 2 (y,T)\{x + y)- 2a dy dx 

K n {x,y)\ci(y,T)\\ci{x,T) - c 2 (x,t)\(x + yY 2a dydx 



o o 

n n—x 




n y 



+ Sn(y)Hx,y)\ci(y,T) - c 2 (y,r)\x a dxdy 







dr 



+ \\u%\\ Y tn 2 ' T Mexp(2tn 2 ' T ML)\\c 1 - c 2 \\ D . 

Since K n (x,y) = for x, y < cr/n, the maximum value that the term (x + y)~ 2a can have is (n/2er) 2<T . Using 
this fact, H5), the definition of M, and Ci,c 2 € D we arrive at 

||G( C1 ) -G(cb) || x _ a „ 

< n 2a Mtexp(2n 2a MLt)\\ Cl ~ c 2 \\ D h~^ +2 ^n 6a a~ 2a M\\ Cl \\\, + Cn e \\c x \\ D \ t 

+ exp (2n 2a MLt) {2~^ 2 ^n 6a (j- 2a M{\\ Cl \\ D + ||c 2 || D )||ci - c 2 \\ D + Cn e \\ Cl - c 2 \\ D ^j t 

+ 1 1 v% (x) 1 1 Y n 2a Mt exp (2n 2<T MLt) \\ci-c 2 \\ D 

< texp (2n 2a MLt) [2 1 - 2a n 6a a- 2a ML 2 t + 2Cn e Lt + 2 1 " 2 " n 6a cj- 2a M L + Cn e + n 2 ° 'ML] \\a - c 2 \\ D 
= texp [2n 2cr MLt) [2 1 ~ 2a n 617 a~ 2a M L(Lt + 1) + Cn 6 (2Lt + 1) + n 2a ML] ||ci - c 2 \\ D , 

from where we can conclude that 

||G(c 1 )-G(c 2 )|| D < 7 ||c 1 -c 2 || fl , 

where 7 = texp (2n 2a MLt) [2 1 - 2a n 6a a^ 2a ML(Lt + 1) + Cn e (2Lt + 1) + n 2a ML] < 1, which completes the 
proof of the lemma. □ 
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Theorem 2.4 Suppose that HI )-H5) hold and u E Y + . Then for each n — 2,3,4, .. . the problem (13)- (16) 
has a unique solution u n with u n (x, t) > for a.e. x £ [0, n] and t £ [0, oo[. Moreover, for all t € [0, oo[ 



xu n (x,t)dx — / xu n (x, 0) dx. 
i) Jo 



(29) 



Proof. From Lemmas |2.1[ |2.3| and the Contraction Mapping Theorem, it follows that there exists a unique 



solution u n (x, t) to ( 16 )-( 18 1 in [0, to}. We proceed now to check that those solutions are non-negative. If we set 

cq=Uq and a = G(c,_i), 
for i = 1, 2, 3, . . ., we find that fixed point iteration gives 

Cj — > u n in L 1 (]0,oo[) as i — > oo, 



and w™ is constructed by positivity preserving iterations, using G given in (17 1 



Let us check now that the mass conservation property (29 1 holds. Multiplying (13) by x and integrating 
with respect to a; on [0, n} we have by |3| and changes of variables and order of integration as in (25 1 



dt 



xu n (x, t) dx 



n n—x 



xK n (x - y,y)u n (x - y,t)u n (y,t) dydx 



xK n (x, y)u n (x, t)u n (y, t) dy dx 



o o 

n n 



xb(x,y)S n (y)u n (y,t) dy dx — / xS n (x)u n (x 7 t) dx 



o x 

n n—x 



(x + y)K n (x,y)u n (x,t)u n (y,t)dydx- / / xK n (x,y)u n (x,t)u n (y,t) dy dx 



o o 

n v 







xb(x,y)S n (y)u n (y,t) dx dy — / xS n (x)u n (x,t) dx 



oo o 

n n—x n n—x 

xK n (x,y)u n (x,t)u n (y,t) dy dx — J J xK n (x,y)u n (x,t)u n (y,t) dy dx 
oo oo 

n n 

+ J yS n (y)u n (y,t) dy - J xS n (x)u n (x,t) dx = 0, 
o o 

from where we have 



— / xu n (xA)dx = 
dt. K ' 



xu n (x 1 t)dx = J xu l (x)dx. 

ii o 
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Now we show that our solution for t £ [0, to] extends to arbitrarily large times, changing variable as we did in 
(25 1 we proceed to obtain a uniform bound 



u n (x,t)x ° dx 



K n (x ~ y,y)u n (x - y,T)u n (y,T)x a dydx 



o L o o 

n n—x 



K n (x,y)u n (x,T)u n (y,T)x a dydx + J J S n (y)b(x,y)u n (y,T)x a dydx 

Ox 



S n (x)u n (x, t)x ° dx 



dr + J Uq(x)x ° dx 
o 



K n (x,y)u n {x,T)u n (y,T)(x + y) a dydx 



o o 

n n—x 



n y 

K n (x,y)u n (x,T)u n (y,T)x~' T dydx + J J S n (y)b(x,y)u n (y,T)x^ (T dxdy 
oo oo 



S n (x)u n (x, t)x "dx 



dr + J Uq(x)x 'dx. 
b 



u n (x,t)x- a dx < 



1 



dr + Uq(x)x a dx 



Making use of the inequality ( 11 1 and the symmetry of K(x, y) results in 

n—a/n n—x 

K(x, y)u n (x, r)u n (y, T)(x~ a + y~ a ) dydx 

cr/n cr/n 
n—a/n n—x 

K(x, y)u n {x 1 T)u n (y, t)(:e~ ct + y^ a ) dy dx 

cr/n a / 7i 

n y n 

+ J J S n {y)b{x,y)u n (y 1 T)x~' 7 dxdy 

cr/n cr/n 
t n y n 

< J J J S n (y)b(x,y)u n (y 1 T)x~ a dxdydT + J u^(x)x^ a dx 

a j n cr/n 
t n y 

< Mn a (j- a J J J u n {y,T)dxdydT+\\v%\\ Y 


t n 

= Mrf<j-° J J yu n (y,T)dydT+\\uZ\\ Y 
o o 

< (n a a- a MT+l)\\u r >\\ Y = L. 



(30) 
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Now we can extend the interval [0,to] to [0, oo[ to complete the proof of Theorem 2.4 By considering the 
operator 



G x (c) {x,t) = - / exp ( - [P x {x, t,c)-Px {x, r, c)] ) 



+ / S n (y)b(x,y)u n (y,t)dy 



K n {x - y, y)c(x - y, r)c(y, r) dy 



dr + u n (x, to) exp (Pi(x, i, c)J , 



with 



Pt(a;,t,c) 



JsT„(a;, y)c(af, r)c(y, r) + S„(x) 



dr, 



we can repeat the above argument to show, that there is a unique non-negative solution u n on [to,ti] where 
tx = 2to. We can extend the unique solution to [0,tj] j = 1,2,3, ... , repeating this process by considering the 
operators 



G j+1 {c)(x,t) = i / exp(- [P j+1 (x,t,c) - P j+1 (x,T,c)]) 



K n (x - y, y)c(x - y, r)c(y, t) dy 



S n (y)b(x,y)u n (y,t)dy 



dr + u n (x, tj) exp (P J+ i(a;, t,c)), 



with 



P j+1 (x,t,c) = 



K n (x, y)c(x, t)c(jj, t) dx + S n (x) 



dr. 



In that way we extend the solution to all of [0, oo[. The argument used to get (29 1 for [0, to] shows that (29 1 
holds for [0, oo[ and thus we have completed the proof of Theorem 2.4 by the arbitrariness of n. □ 



2.2 Properties of the solutions of the truncated problem 



Lemma 2.5 Let u n a solution of the truncated problem { 13)- ftl4\ l- Then for < a < 1 and n = 1,2,... we 
obtain the inequality 



dt 



u n {x,t)x- a dx< / S n (y)b(x,y)u n (y,t)x- a dy dx. 



x 



Proof. Multiplying equation (13 1 by a; a and integrating w.r.t x from to n we have 

n n x 

d 



. u n (x, t)x a dx = - 
dt K ' 2 



o o 

n n—x 



K n {x - y, y)u n (x - y, t)u n (y, t)x a dy dx 
K n {x, y)u n {x, t)u n (y, t)x~ a dy dx 



o o 

n n 



+ J J S n (y)b(x,y)u n (y,t)x~ a dydx- J S n (x)u n (x,t)x~ a dx. 

Ox 
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Changing variables as we did in ( 25 1 we get 

n n n — x 

~ f u n (x,t)x~ a dx = - [ [ K n (x,y)u n (x,t)u n (y,t)(x + y)- a dydx 



dt y ' ' 2 



o o 

n n—x 



K n (x,y)u n (x,t)u n (y,t)x a dydx 



o o 

n n 



+ / / S n (y)b(x,y)u n (y,t)x- a dydx- / S n (x)u n (x,t) X - a dx 



Now by using inequality (11) together with the definition and symmetry of K n (x, y) 

n n n—x 

j t [ u n (x,t)x- a dx < 1 2 -( Q+1 ) f ( K n (x,y)u n (x 7 t)u n (y,t)(x- a + y- a )dydx 



o o 

n n— x 



-\j J K n (x,y)u n (x,t)u n (y,t)(x- a +y- a )dydx 



n n n 

+ J J S n (y)b(x,y)u n (y,t)x- a dydx- J S n (x)u n (x,t)x- a dx 

Ox 
n n—x 

= \ (2-(« +1 ) - l) J J K n (x,y)u n (x,t)u n (y,t)(x- a +y- a )dydx 
o o 

n n n 

+ J J S n (y)b(x,y)u n (y,t)x~ a dy dx - J S n (x)u n (x,t)x~ a dx. 

Ox 

Now we can eliminate the negative terms and obtain 

n n n 

j t Ju n (x,t)x~ a dx< J I S n (y)b(x,y)u n (y,t)x' a dydx, 

Ox 

which complete the proof of the theorem. □ 
In the rest of the paper we consider for each u n their zero extension on K, i.e. 



u n (x,t) 



u n (x,t) 0<x<n, te [0,T], 
x < or x > n. 



For clarity we drop the notation • for the remainder of the paper. 

Lemma 2.6 Assume that Hl)-H5) hold. We take u n to be the non-negative zero extension of the solution to 
the truncated problem found in Theorem\2.J\ Fix T > and let us define 



L(T) = (e NT (N + 1) + e CT (C + 1) + l) Klb 

Then the following are true: 

i) We have the bound 

(l + x + x- 2a )u n (x,t)dx < L(T) for all t e [0,T]. 

ii) Given e > there exists an R > 1 such that for all t € [0, T] 

sup I f(l + x-" 7 )^^, t)dx\ < e. 
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iii) Given e > there exists a S > such that for all n = 2, 3, . . . and t G [0, T] 
\l + X -^( X ,t) dx <e W hene V er »(A)< 5. 



Proof. Property i) By Lemma 2.5 for a — we have 



u n (x,t)dx< / S n (y)b(x,y)u n (y,t)dydx. 



dt 

Ox 

Changing the order of integration in the right hand side and making use of ([4} and H4) , we get 

n n y 

d 



dt 



u n (x,t)dx< I S n (y)u n (y,t) / b(x,y)dxdy 



N J y°u n (y,t) dy + Nj y»u n (y,t)dy. 
o 1 



As 9 is considered to be in [0, 1[, by using the mass conservation property (29 1 we have the estimate 

n n 

d 



dt 



II: 



(x, t)dx<N u n (y, t) dy + iV||M ||y. 





Integrating respect to time it becomes 

n t n n 

u n (x,t) dx < N J J u" (y,T)dydT + N\\u a \\ Y - J u£(x)dx 



t n 



<NJ J u n (y , t) dy dr + (N + 1)\\u \\y , 
o o 

from which we obtain by the Gronwall's inequality, see e.g. Walter [Til page 361], for P(t) = /„" u n (x, 
u n {x)dx<e m {N + l)\\u \\ Y , f G [0,T]. 



Computing now the term with the weight x using Lemma 2.5 for a = la we have 



~J u n (x,t)x- AcT dx< J S(y)b(x,y)u"(y,t)x-^dydx. 

a: 

Changing the order of integration and using H5) we get 

n n y 

-2a, 



j I u n {x,t)x- 2,y dx 



S{y)b{x,y)u n (y,t)x- 2a dxdy 



o 







< Cj y e ~ 2 °u n { Vl t)dy 
o 

1 n 
<cj y-^u n (y,t)dy + C J yu n {y,t)dy 

1 

71 

<C J x- 2a u n (x,t)dx + C\\u \\ Y . 
b 
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From this inequality we find as above that 

u n (x,t)x~ 2rT dx < e ct (C + l)\\u \\ Y , tG [0,T]. 



Now, by the mass conservation property (29 1, by (31 1, and (32) we obtain 

n n n 

,i + x + *-)„•<,. ., * - / «•■(,, o * + / „■(,. „* + / o * 



< e Nt (N + l)\\u Q \\ Y + ho||y + e ct (C + 1)\\u \\y 

< (e NT (N + 1) + e CT (C + 1) + 1) hollr =: L{T). 

Property ii) Choose e > and let R > 1 be such that R > 2 ^ M °H y . 



Then using ( 29 1 we 



(1 + x^ a )u n (x,t)dx = / -M n (x,t)rfx+ \ x~ a - I -^u n (x,t)dx 



< — / xu n (x,t)dx + 



1 



< 



1 1 



sit™ (a;, i) ciir 



xu n (x, t) dx 



< — J xu n (x, t) dx 



2 2 
< ^||«oI|y < ^ll«o||y < c 



Property iii) By property (m) we can choose r > 1 such that for all n and i € [0, T] 



(1 + x -<7 )u n (a:,i) dx < 



By using H5) and ([4} we find that 
v i 
b(x,y)x~ a dx = / y)x~ a dx + I b(x,y)x~ a dx 



< j b{x,y)x a dx + j b(x,y)dx 
o 1 

< J b(x,y)x- 2a dx + J b(x,y)dx <Cy- 2a + N. 
o o 

From this we get for y G [0, r] 



S(y)b(x 1 y)y a x ° dx < j y b(x,y)y a x dx 
o 



< r e+2a / b{x,y)x- a dx < Cr e + Nr e+2 ° 



Let xa denote the characteristic function of a set A, i.e. 
, s f 1 if x e A 
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and set 

«(r) = - max (1 + y CT )(l + x)\\ + y) A . (35) 

Z 0<a:<r 
0<y<r 

Then, by the absolute continuity of the Lebesgue integral there exists a 81 > such that 

y 

sup xa(x + z)S{y)b{x 7 y)y 2a x- a dx < ——— -, (36) 

o< z <r J 8L(T)T exp ln(r)L(T)T) 

0<y<r 

whenever ^4 c]0,oo[ with /j.(A) < Si. Also by ^ and the abolute continuity of the Lebesgue integral there 
exists a S 2 > such that 

SU P / XAnlO r]( x + -zWx, w) dx < — - — -, , - — - — ; — r , (37) 

0<z<rJ 1,1 ~ 8L(T)Texp(K(r)L(T)T) , y ' 

o 

for fi(A) < 82- Now, ||uo||y < L leads to 

00 00 1 00 

J(l + x-nu (x )d x < Ju Q (x )d x + Jx-°u Q (x )d x + Jx-°u Q (x) d x 

1 

00 1 



< Ju a (x)dx + Jx^u (x)dx + Ju (x)dx 

1 

< 2||w ||y < 2L. 



Let us define for all n = 1,2,3,... and t e [0, T] 

00 

f n (A,t) = sup / XAn[a.r]{x + z){l + x- a )u n {x,t) d x. 

0<z<rJ 


By the absolute continuity of the Lebesgue integral there exists a £3 > such that 

00 

/"(AO) = sup ( XAn%r}{x + z){l + x-°)u n {x) d x< - \ (38) 

o< z < r J 4exp (k{r)L(T)T) 

for n(A) < 83. Now we take S = min(5i, 82, 83) and multiply ( 13 1 by (1 + x~ a )xAn[Q.r] (% + z ) f° r z € [0, r]. This 
we integrate from to t w.r.t. s and over [0, oo[ w.r.t. x. Using the non-negativity of each u n and fi(A) < 8 we 
obtain 



00 

J XAn[o,r] {x + z)(l + x~ a )u n (x, t) dx 
'0 

t 00 00 

-2/// XAn[o,»i(z + ^)X[o,*]n[o,r](y)(l + ^~ CT )^n(^-y,^ 
000 

i 00 n 

+ y J XAn[o,r](x + z) J S n (y)b(x,y)(l + x~ a )u n (y, s) dydx ds 

00 k 

00 

+ y XAn[o,r](x + z)(l + x~ a )uo(x)dx. 
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By changing variable as we did in ( 25 1 we get 



1 

< 

~ 2 



XAn[o,r]{x + z)(l + x a )u n (x,t)dx 

XAn[0,r]( x + V + z )X[o,x+y]n[o,r](y)[^ + + yr a ]K n (x,y)u n (x, s)u n (y, s) dydxds 



o 

t oo oo 





t oo 



XAn[o,r# + z) / S n (y)b(x,y)(l + x a )u n (y, s) dy dx ds + f n (A,0) 







l2i+l22 + f n (A,0). (39) 



Now, we estimate the integral terms I21 and I22 in (39 1. Working with I21, by using the estimation of K(x,y), 



then taking 1 + (x + y) a < 1 + y a and x a < 1 + x " we have 

t 00 00 

hi<^J J J XAn[o.r]{x + y + z)x[a.x+ y \n[o,r]{y) [l + {x + y)~ a ] (± + x) x (1 + y) x (xyY 
000 

■«™ (x, s)u n (y, s) da; dy ds 

t 00 00 

1 f f f , s v A/-, „.\A 



< 

~ 2 



000 



< 
~ 2 



i OO OO 

1 



000 



XAn[o,r](a; + y + «)X[o^+«]n[o,r](y)(l + y a )(l + x) x (l + y) x (l + x a )y-° 
■u n (x, s)u n (y, s) dx dy ds 

XAn[o,r](x + y + «)X[o J «+»]n[o,r](y)(l + + ^) A (1 + 2/) A (l + ^ a )y~ 2 " 
•it™(x, s)u n (y, s) dx dy ds. 



By using the definition (35 1 of n(r) and property (i) we obtain the following estimation for I21 



hi < k(r) J J u n (y,s)y J XAn[o,r](x + V + z)(l + x <T )u n (x,s)dxdyds 
00 

t r oo 

< k(r) u n (y,s)y~ 2 ' 7 sup / XAn[o.r] (x + w)(l + x~ a )u n {x, s) dx dy ds 

J J 0<ui<rJ 


t 

<k(r)L(T)Jf n (A,s)ds. (40) 



Working now with the integral term I22, by changing the order of integration we get 

t r n 





t r y 



I22 = I I XAn[o,r](x + z) J S n (y)b(x,y)(l+x a )u n (y, s) dy dx ds 

XAn[o.r] (x + z)S n (y)b(x, y)(l + x _<T )u"(y, s) dx dy ds. (41) 



000 

t n r 
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As r < y in the second integral term of (41 1 we find that 

try 

hi< J J J XAn[o,r](s + z)S n (y)b(x,y){l + x~ a )u n (y,s) dxdyds 



t n y 



r 
t n y 




t n y 





t n J/ 



XAn[0T](i + z)S n (y)b(x,y)(l + x _CT )u"(y, s) dxdyds 
XAn[o,r] {x + z)S n (y)b(x, y)(l + x~ a )u n (y, s) dx dy ds 
XAn[o,r](x + z)S n (y)b(x,y)u n (y, s) dxdyds 

XAn[o,r](x + z)S n {y)b{x,y)x~ a u n {y,s) dx dy ds. 



ooo 



By using (37 1 and property (i), we estimate the first integral term on the right hand side of (42 1 by 

t n y 

XAn[o,r] + z)b(x, y)S n (y)u n (y, s) dx dy ds 



ooo 

t n 



II 



< / / S n (y)u n (y,s) sup / XAn[o,r](x + z)b(x,y) dx dy ds 

J J 0<z<rJ 




< 



< 



8L(T)Texp( K (r)L(T)T) 



y°u n {y,s)dy+ / y»u n (y,s)dy 



ds 



8L(T)T exp ( K (r)L(T)T) 



ds < 



8exp( K (r)L(T)T)' 



u n (y,s)dy + J yu n (y,s)dy 

.0 1 

Now, using (36 1 and (i) we also estimate the second integral term on the right hand side of ( [42] ) by 




t n 



y 



< 



u n (y,s)y a sup XAn\o.r]{x + z)S n (y)b(x,y)x a y a dxdyds 

0<Z<7- J 

0<y<r 



< 



&L(T)T exp U{r)L{T)T) 



u n {y,s)y- 2a dyds< 



8exp(n{r)L(T)T) 



This brings together with the estimation (43 1 of the first integral term of (|42j) 
hi < 



4exp(«;(r)L(r)T) 



Now we insert into a39\ the estimates ( 40 1 , ( 44 ) . and ( 38 1 to arrive at 



XAn[o,r]{x + z){l + x-' J )u n {x,t)dx < k{r)L{T) / f n (A lS )ds + 



2exp(k(r)L(T)T) ' 

u u 

Since the right hand side is independent of z we may take sup 0< ^ <r on the left hand side to obtain 



P(A,t)<k(r)L(T) / P(A,s)ds + 



2exp(k(r)L(T)T)' 
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By using Gronwall's inequality see e.g. Walter [l4j page 361], we get 
f n (A,t) < eexp (k(r)L(T)T)/ (2exp (k(r)L(T)T)) = e/2. 



(45) 



Thus, by p§) and (45 1 for every n = 1, 2, 3, . . . and £ e [0, T] 



(1 + ac a )u n (x,t)dx = j XAn[o,r](l + a; ct )u™(:e, i)dx + / XAn[r,oo] (1 + a; a )u n (x,t)dx 

AO 

oo 

< / (l + a;- CT )u"(a;,t)da;<e/2 + e/2 = e 



whenever n(A) < S. 

This completes the proof of Lemma |2.6| □ 
Let us define v n (x, t) — x~ a u n (x, t). Due to the Lemma 2.6 above and the Dunford-Pettis Theorem [3] page 

274], we can conclude that for each t £ [0, T] the sequences (M n (i)) ngN and ( w "(*)) neN are weakly relatively 

compact in L 1 (]0, oo[). 



2.3 Equicontinuity in time 

Lemma 2.7 Assume that Hypotheses 1.1 ftoW. Tafce (u n ) now £0 &e £/ie sequence of extended solutions to the 
truncated problems (13)-u6V found in Theorem 2.4 and v n (x,t) — x~"u n [x,t). Then there exists a subsequences 
(it nfc (t)) and (i> nj (t)) 0/ (u n (i)) N and ngN respectively such that 

u 7lk (t)^u(t) in L 1 (]0,oo[) as n* 00 
i; ni (i)^n(t) m L 1 (]0,oo[) as n t — s- 00 

uniformly for t £ [0,T]. XTiis convergence is uniform for all t £ [0;T] giving u, v € Cb([0, T]; L 1 (]0, oo[)) = 
{77 : [0,oo[— >• L 1 (]0, oo[), 77 continuous andi](t) bounded for allt > 0}. 



Proof: Choose e > and 6 L°°(]0, oo[). Let s,t G [0, T] and assume that t > s. Choose a > 1 such that 

||<PlU°°(]o,<x>[) < e/2. 



(46) 



Using Lemma |2.6[i), for each n, we have 



|M™(a;,i)-M n (x,s)|dx< - / x \u n (x,t) + u n (x, s)\ dx < 2L(T)/a. 

a 



(47) 



By multiplying (13 1 by cf> and integrating w.r.t. x from to a as well as from a to 00, w.r.t. r form s to t and 



using (46 1, (47 1 and t> s we get 



< 



< 



4>(x) [u n (x, t) - u n (x, a)] dx 



(x) [u n (x, t) - u n (x, a)] dx 



(x) [u n (x, t) - u n (x, a)] dx 



H0l|i°°C]O,co[) 



K„(x - y,y)u n (x - y , t)u 71 (y , t) dy dx + / K n (x,y)u n (x,T)u n (y,T) dy dx 











+ j J b(x,y)S n (y)u n (y,T)dydx + j S n (x)u n (x,T) dx 

Ox 



dr + e/2. 
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Changing variables as we did in ( 25 1 we have 



<j)(x) [u n (x, t) - u n (x, s)] dx 



||0||l°°(]o,ooD 



K n (x,y)u n (x,T)u n (y,T)dydx+J J K n (x,y)u n (x,r)u n (y,T) dy dx 
oo oo 



+ J J b(x,y)S n (y)u n (y,T)dydx + j S n {x)u n (x,T) dx 

Ox 



dr + e/2. 



Using the definition (12 1 of K n (x,y) and S n (y) and the estimation of K(x,y) we obtain 
(x) [u n (x, t) — u n (x, s)] dx 



< ||^l|i~(]0,oo[) 



a—cr/n a—x 



(1 + x)\l + y) x {xy)- a u n {x, r)u n (y, r) dydx 



a/n a/71 

a n—x 



+ I I {l + x) x (l + y) x {xy)- a u n {x,T)u n {y,T)dydx 

a/n a/n 



b(x, y)S(y)u n (y, r) dydx + J x u n (x, r) dx 

a/n x a/n 

IMU-(]o,«,D / (Mir) + I 3 2(r) + I 33 (t) + I 3i (T)) dr + e/2. 



dr + e/2 



(48) 



Now, we estimate the terms I 3 i(t), I 3 2{t), ^33(1"), and I 3 i(r) in (48 1. By using Lemma 2.6 i) the first term can 
be estimated by 



hi(r) <-(! + «> 



a-a/n a —x 



2 A 



1 



{xy)-"u n {x, T)u n (y, r) dydx<-{\ + a) 2X L(T) 2 . 



a/n a/n 

In order to estimate the second term, we define 

,1 if < A < 1 
2 A -! if A>1. 



(49) 



(50) 



Then, by using inequalities pi and (10 I for p = A and Lemma [2. 6| (i) we find that 



/32(r)<Ci J I (l + x x )(l + y x )(xy)- a u n (x,r)u n (y,T)dydx 

a/n a/n 

a n—x 

<d(l + a A ) J J {l + y x ){xy)-°u n {x,T)u n {y,T)dydx 

cr/n ajn 
a n—x 

<Ci(l + a A ) J J (y- a +y x -' 7 )x- a u n {x,T)u n (y,T)dydx<2C 1 (l + a x )L(T) 2 

a/n a/n 



(51) 
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We also find by using H4). Lemma 2.6 (i), and Q that 



I 33 (t)= / J b(x,y)S(y)u n (y,T)dydx 

a/n X 

ay n a 

b(x,y)y e u n (y,T)dxdy + J J b{x 1 y)y e u n {y,r)dxdy 

a/n a fn & a/n 



<N j y a u n (y,r)dy<NL(T). 

a/n 

Now, by using Lemma |2.6[i) we have 



(52) 



I 34 (t)= / x u n (x,r)dx < L(T), 

cr/n 



(53) 



which together with (49l-(52l brings (48) to 



f>(x) [u n (x, t) - u 1l (x, a)] dx 
Q(l + a) 2A + 2d (1 + a x )^j L(T) 2 + (N + l)L(T) 



(t-a)H^IU- +e/2 < e, 



(54) 



whenever (t — s) < S for some 5 > sufficiently small. The argument given above similarly holds for s < t. 
Hence ( 54 1 holds for all n and \t — s| < S. Then the sequence (u n (t)^ nefsj is time equicontinuous in L 1 (]0, oo[ 



Thus, (u n (t)) lies in a relatively compact subset of a gauge space f^. The gauge space is i 1 (]0, oo[) equipped 
with the weak topology. For details about the gauge space, see Ash [l] page 226]. Then, we may apply a version 
of the Arzela-Ascoli Theorem, see Ash fll page 228], to conclude that there exists a subsequence (u nk ) k(£N such 
that 

u nk (t)^-u(t) in fix as n k -> oo, 

uniformly for f G [0,T] for some u G C([0, T]; L x (]0, oo[)). 

Now let us consider v n {x, t) — x~ a u n (x, t) where we have to deal with a stronger singularity at 0. 
We take e > 0, <j>, s and t as they were defined before. Using Lemma [2T6| for each n, we get using a > 1 chosen 



to satisfy ( 46 ) 



\v n (x,t) — v n (x,s)\dx = / \x °u n (x,t) — x a u n (x, s) I dx 



< - I x l - a \u n (x,t)+u n (x,s)\dx 



< - / x\u n (x,t) +u n (x,s)\dx < 2L(T)/a. 
a 



(55) 
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By using (131, (46 1, (55 1, for t > s and the definition of v n {x) we obtain 



0(x) [v n (x,t)-v n (x,s)} dx 



< \<f>{x)\[v n (x,t)+v n (x,s)]dx + e/2 



< ||0|U-(]O,oo[) 



K n (x - y,y)u n (x - y, r)u n (y, t)x a dydx 



o o 

a n—x 



+ J J K n (x,y)u n (x,T)u n (y,T)x a dydx 



an a 

+ \l b(x,y)S n (y)u n (y,T)dydx+ / S n (x)u n (x, r) dx 



x 







= IMU-ao,ooD / ( J 4i(r)+l42(r) + 743(T)+l44(r))dT + e/2. 

s 

A change of variables in the first integral gives 



dr + e/2 



1 



2f„0, y)u n (x, r)u n (y, t)(x + y) a dy dx 



o o 



Taking y = in the term (x + y) " we find that 



hi(r)<~ / K n (x,y)u n (x,T)u n (y,T)x a dydx 



o o 



Working in a similar way as we did in (49 1 and (51 1 we find the estimations 

a—cr/n a—x 



hi(r) 



1 



(1 + x) A {l + y) x x- 2a y- a u n {x, T )u n (x, r) dy dx 



a I n a/r 



a— a jn a -x 



x- Za y-°u n {x 1 T)u n {x,T)dydx < -(1 + a y x L{T) 



a/n a/n 



and 



a n — x 



Ia2{t)<C 1 J I (l + x A )(l+y A )x-^y- a u n (x,T)u n {y,T)dydx 

a/n a/n 



a n—x 

<Ci(l + a A ) J J (y- a + y x -' 7 )x- 2 ' T u n (x,T)u n (y,T)dydx<2C 1 (l + a x )L(T) 2 

a/n (7 jn 

Now changing the order of integration in I^{r) we have 

a n 

7 43 (r)= j j 'b(x 7 y)S(y)u n ( yi T)x-°dydx 

a/n x 

ay n a 

b(x, y)y 9 u n (y, r)x~ a dx dy + J J b(x, y)y e u n (y, r)x~° dx dy. 

a jn a/n & a jn 
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By using ( 34 1 and H4) we get 



hs(T)< J (Cy~ 2a + N)y»u n ( yi T)dy + J {Cy- 2a + N)y H u n (y,T) dy 

cr/n a 
(Cy- 2 ° + N)y e u n (y,T)dy 

cr/n 

n n 

-C J y e - 2 °u n (y,T)dy + N J y e u n (y,r) dy < (C + N)L(T). 

cr/n cr/n 



Using H4) and Lemma 2.6 (i) we obtain 



hi{r)= J x u n (x, r)dx < L(T) 

cr/n 



(59) 



(60) 



which together with (56l-(59l gives the estimation 



x) [v n (x, t) — v n (x, s)] dx 

1(1 + a) 2X + 2Ci (1 + a x ) ) L(T) 2 + (C + N + l)L{T) 



(i-s)||0|| L ^ (] o,oo[)+e/2. 



We can use now the same argument used for u n to conclude that there exists a subsequence (i | " fc ) fe6N such that 

v Uk (t) — > v(t) in Q as n k -> oo, 

uniformly for t G [0,T] for some v € C([0, T]; L^QO, oo[)). 

Since T > is arbitrary we obtain u,v g Cs([0, oo[; L 1 (]0, oo[)). 

Lemma 2.8 For v n (-,t) defined as before, we have 

V n (-, t) — v(-, t) where v(x, t) = x- a u(x,t) for all t € [0,T] in L^QO,^). 



□ 



Proof. By Lemma 2.7 we know that v n (t) v(t) in L 1 (]0,oo[) as n — > oo uniformly for < € [0, T]. 
Then, we just need to prove that £) = x _<T u(a;, i) 
By definition of weak convergence we have 

<p(x) [v n {x,t) -v(x,t)]dx-> for all ^ G L°°(]0, a]). 
As x CT e L°°(]0,a]) we obtain 

a 

^(x) [^"(x, t) - x a f (x, *)] dx = J <p(x) [u n (x, t) - x a v(x, t)] dx -> for all 93 e L°°(]0, a]). 

Since u" — 1 u we have due to the uniqueness of the limit of weak convergence, v(x, t) — x~ a u(x, t). □ 

3 Existence Theorem 

3.1 Convergence of the integrals 

In order to show that the limit function which we obtained above is indeed a solution to ([TJ-Q, we define the 
operators Mf , M h i = 1, 2, 3,4 



Mf(u")(x) = lfK n {x-y,y)u n (x-y)u n {y)dy M^x) = \] K(x - y,y)u(x - y)u{y)dy 



n—x oo 

M 2 n (0(x)= / K n (x,y)u n (x)u n (y)dy M 2 (u) (x) = J K(x, y)u[x)u(y) dy 



n oo 

M$(u n ){x) = Jb(x,y)S n (y)u n (y)dy M 3 (u)(x) = J b(x , y) S (y)u(y) dy 



M%{u n )(x) =S n {x)u n {x) 



Mi{u){x) = S(x)u(x), 
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where u G i 1 (]0, cxd[), at e]0,oo[andn = 1,2,.... Set M" = M™ - Af™ + M 3 " - and M = M X -M 2 +M^-M A . 

Lemma 3.1 Suppose that (u n )„ eN C Y + , u G ^ + where ||u"||y < i, < Q, u n ^ u and v n — 1 u in 
L 1 (]0,cx)[) as n — > oo, /or u = x~ a u and v n — x~ a tt n . Then for each a > 

M"(u n ) M(u) in I^QO.aQ as n -> oo. 



Proof: Choose a > and let G L°°(]0,oo[). We show that Af"(w") M;(u) in ,L 1 (]0, a[) as n -> oo for 
i = 1,2,3,4. 

Case i = 1: For u G Y + and a; G [0, a] we define the operator g by 



= - / <fi(x + y)K(x,y)(xy) a v(y) dy where v = x a u. 
o 

For a.e. x G [0, a] the function defined by 

<Px{y) ■■= ^X[o, a -x](y)4>(x + v)K{x,y){xyY < ^X[a,a-x](y)<t>{x + y){l + x) x (l + y) x , 

where x denotes the characteristic function and the estimate is due to H3), is in L°°Q0, oo[). Since v n — 1 v in 
L^OjOoI), it follows that 

g(v n )(x) -> g(v)(x) for a.e. xe[0,a]. (61) 

OO 

Also, as ||w n ||y < L implies J v n (x) dx < L we have 





a— x 



< 



(x + y)(l + x) A (l + y)V(y)d 2 / 



This holds analogously for g(v). Thus, both, g(v n ) and g(v) are in L°°([0, a]) with bound 

llff(«")IU-([0,a]) + lls(«)IU~([0,a]) < + a) 2Q H0llL~([O,a])(£ + Q)- 



(62) 
(63) 



It follows by (61 1 and Egoroff's Theorem |TJ page 94] that 

g(v n ) — > g(v) as n — > oo almost uniformly in [0, a]. (64) 

Remember that almost uniform convergence means that for any given 5 there exists a set EC [0, a] such that 
fj.(E) < S and g(v n ) — > g(v) uniformly on [0, a] \ E as n — ¥ oo. 

By Lemma 2.6 [Hi), since v n — 1 u in L 1 (]0,oo[), there is a 5 > such that for all n 



u"(a;) < e/ [(1 + a) 2a ||0|| i =o([ OiQ ] ) (L + Q)] whenever fi(E) < 8. 



(65) 



By the Holder inequality we obtain using (|63| and (|65|) 

[ 5 (*;™)(z)-^)(x)K(x)dz 



[ff(« B )(x)- 5 («)(x)] 



'[0,a]\B 

< llffK) -5(«)|U«>([0,o]\B) 

< ||fl(« n ) -5(u)||L»([0,a]\B) 



[g(v n )(x)-g(v)(x)]v n (x)dx 



v n (x)dx+[\\g(v n )\\ L ~ (E) + \\g(v)\\ L ~( E) } / i^dz 

[0,a]\B JE 



[0,a]\E 



v n (x) dx + — < e for n > uq. 



2o 



Since e > was arbitrarily chosen the almost uniform convergence of g(v n ) leads to 
[g(v n )(x) - g(v)(x)]v n (x) dx as rn-oo. 



(66) 



Also, since g(v) £ L°°([0, a]) is bounded independently of n by ( 62 1 and v n — 1 v in L 1 (]0,cxd[) as n — > oo 



g r (u)(a;) [w n (a;) — v(x)) dx 



— > as n — > oo. 



(67) 



Now, since g(v n ) € L°°([0, a]) and d" e ^ 1 ([0, oo[), by the absolute continuity of the Lebesgue integral, we have 



< ||s(« n )IU~ a0lO]) 



a — a/r 



as 7i — > oo, 



g{v n ){x)v n {x) dx 

a — a/n 

In the same way we get 

J g(v n )(x)v n (x)dx 
o 

and 

a—a/n a/n 

J J (j>(x + y)K(x, y)u{x)u{y) dy dx 

a/n 

a—a/n a/n 

J J cj){x + y){l + x) x {l + y) x v{x)v{y)dydx 

a/n 

Now, we show that 



v n (x) dx 



— > as n — > oo, 



< 



— » as n — >• oo. 



(68) 



(69) 



cj)(x)M?(u n )(x)dx 



i x )^ I K n{x - y,y)u n (x - y)u n {y)dydx 



o o 

a a—x 





a—a/n a -x 



1 

2 

1 
2 

1 

2 

<r/n 
a—a/n 



<j>(x + y)K n {x 1 y)u n (x)u n {y) dy dx 



4>(x + y)K(x, y)u n (x)u n (y) dydx 



a/n a/n 
a—a/n a—x 



a—a/n a jr. 



(x + y)K{x,y)u n {x)u n {y) dy dx - - 



<j>(x + y)K(x, y)u n (x)u n {y) dydx 



a/n 



1 



a— a /n a/n 



g(v n )(x)v n (x)dx- | / / <j>(x + y)K(x,y)u n (x)u n (y)dydx. 

a/n a jn 

In a similar way we also find that 

a a 

x)Mi(u)(x) dx — J g(v)(x)x~ a u(x) dx — J g(v)(x)v(x) dx. 





(70) 



(71) 



2(> 



Now. it follows from (70 1 and (71) that 



(f>(x) [M?(u n )(x) - Mi(u)(a;)] dx 



a— a jr. 



g(v n )(x)v n (x) dx 



1 



a— a I n cr/r< 



(x + y)K(x,y)u n (x)u n (y) dy dx — / g(v)(x)v(x) dx 



a I n a/n 

a 

g(v n )(x)v n (x) dx — I g{v){x)v(x) dx 



g(v n )(x)v n {x)dx 



a— cr/r 

a— a I n a jn 



g(v n ){x)v n (x) dx- 



1 



<j)(x + y)K(x, y)u n (x)u n (y) dy dx 



< 



a/a 

g(v n )(x)v n (x) dx — I g(v)(x)v(x) dx 



+ 



g(v n )(x)v n (x) dx 



a— a I r. 



tr/r 



g(v n )(x)v n (x)dx 



a—a/n a jr. 



(x + y)K(x, y)u n (x)u n (y) dy dx 



a jn 



By addition and subtraction of the term J g(y){x)v n (x) dx in the first term of the above inequality, it results 

o 

that 



(j>{x) [M[ L {u n ){x) - M x {u){x)] dx 



< 



g(v)(x) [v n (x) — v(x)] dx 



[g(v n )(x)-g(v)(x)]v n (x)dx 



g(v n )(x)v n (x)dx 



a—a/n 



g(v n )(x)v n (x) dx 



a—a/n a jn 



4>{x + y)K{x, y)u n (x)u n {y) dy dx 



a/n 



Now, by ( 66 1-( 69 ) and taking n->cowe have 

(j>{x) [M"(u n ){x) - Mi(u)(x)} dx as n oo. 

It follows, since <j> is arbitrary, that 

M?{u n ){x) M x {u){x) in L X (]0,a[) as n -> oo. 



(72) 



(73) 



Case i = 2: For every e > and C\ defined by ( 50 ) we can choose r\ large enough, due to the negative exponents, 
such that for L, Q from our assumptions 



CilMU~([0,«]) [(2ry-( 1+ff ) + 7 ? A — 1 )(i 2 + Q 2 ) 
Redefining the operator g for u G Y + and x G [0, a] by 

g(v)(x)= c/){x)K(x,y)(xy) a v(y)dy. 



e 

<3- 



(74) 
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For a.e. x e [0, a] the function defined by 

<Px{y) ■= ^X[o, v ](y)<P(x + y)K(x,y)(xy) a 

where, as before, x denotes the characteristic function, is in L°°(]0, oof). Using a similar argument as the one 
that was used in ( 61 )-( 66 1 it can be shown that also for the above redefined g (66 1 and (67) hold. By H3) and 
( 50 1 we have 



<p(x)K(x, y) [u n (x)u n (y) — u(x)u(y)] dy dx 



71 



< 



dl I \<f>(x)\ [(xy)-° + x x -°y-° + x- a y x -° + (xy)^] [u n {x)u n {y) + u{x)u{y)\ dydx 
o n 



< Ci||<^||i/x.([o, a ]) 



(xy) a [u n {x)u n {y) + u(x)u(y)} dydx 



.0 V 



a oo 

X X - a y- T \u , '[.r)u n (u) + »(./■)»(//)] dyd.r 

o n 

a oo 

+ J J x~ a y X - a [u n (x)u n (y) + u(x)u(y)} dydx 
o n 

a oo 

+ J J(xy) x -° [u n (x)u n (y)+u(x)u(y)]dydx 

r. 



We can estimate the first integral term of ( 75 1 as follows 



(xy) a [u n (x)u n (y) + u(x)u(y)\ dydx 



?? 

OO 

< 



'/ 



x a [u n (x)u n (y) + u(x)u(y)] dx + x a \u n (x)u n (y) + u(x)u(y)] dx 



dy 



< 



< 



a 

~ x [u n (x)u n (y) + u(x)u(y)} dx + J x [u n (x)u n (y) + u(x)u(y)] 



dx 



dy 



(x 1 + x) [u n (x)u n (y) + u(x)u(y)] dx 



o 

oo 



dy 



< rr (1+ff) / y [Lu n (y) + Qu(y)} dy < ^(^(L 2 + Q 2 ). 
v 

In the similar way we have 

x x - a y- a [u n (x)u n (y) + u(x)u(y)] dydx < r)-^ +a \L 2 + Q 2 ), 



(75) 



(76) 



(77) 



7} 



x-°y x -° [u n {x)u n {y) + u(x)u(y)] dydx < ■q-^ +(J - x \L 2 + Q 2 ), 



r\ 



and 



(xy) X ~ a [u n {x)u n (y) + u(x)u(y)} dydx < i]~( 1+cr ~ x \L 2 + Q 2 ). 



(78) 



(79) 



7] 
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By ( 76 1-( 79 1 and (74 1, (75) becomes 



4>{x)K{x,y) [u n (x)u n (y) — u(x)u(y)] dydx 



o n 



< 



Now, using Lemma 2.6 [i) and the absolute continuity of the Lebesgue integral, we have 

a <r/n 



4>(x)K(x, y)u n (x)u n (y) dydx 







— ¥ - for n larger than some no , 
3 



(80) 



(81) 



and 



a/n n- 



cj)(x)K(x 1 y)u n (x)u n {y) dy dx 



cr/r 



— > — for n > jiq. 



(82) 



Also, proceeding as before, for n > a we have 



(j)(x)K(x,y)u n (x)u n (y) dydx 



< 2 



(L 2 + Q 2 )n\\ L ~ ao , a]) . (83) 



From ( 80 (-( 83 1 together with the analogues of (66 1 and (67 1, for n > a 



<j>(x) [M$(u n )(x) - M 2 (u)(x)] dx 



(j>(x)K(x,y)u n (x)u n (y) dy dx — / / cf>(x)K(x,y)u(x)u(y) dy dx 



o o 

a cr/n 





a/n n - 



4>(x)K(x,y)u n (x)u n (y) dydx - 



4>{x)K{x, y)u n (x)u n (y) dydx 

a/n 

a rj a oo 

j)(x)K(x,y) [u n (x)u n (y) - u(x)u(y)} dy dx + J J <f>(x)K(x,y) [u n (x)u n (y) - u(x)u(y)] dydx 

)) 

a cr/n v/rin-x 

4>{x)K(x,y)u n (x)u n (y)dydx- I I (j)(x)K(x 1 y)u n (x)u n (y) dy dx 



o o 



cr/r 



< 



4>{x)K(x,y)u n (x)u n (y) dydx 



[g(v n )(x)-g(v)(x)]v n (x)dx 



g(v)(x) [v n (x) — v(x)) dx 



+2\{n-a)- ( - 1+ ^ + {n-a)-^ 1+a - x ^ {L 2 + Q 2 )\\<P\\ L ~ {[QM) <e as n -> oo. 
Therefore, since <p and e are arbitrary, we conclude that 

M%(u n ){x) -± M 2 {u){x) in L x Q0,a\) as n -> oo. 



(84) 



Case i = 3: As 9 G [0, 1[ in H4) we can choose r > a such that for N from Q and L, Q from our assumptions 



(85) 
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By a change of the order of integration and (|4| we have 



<j)(x)b(x,y)S(y) [u n {y) - u(y)} dydx 



a/n r 



(/>{x)b(x,y)S(y) [u n {y) - u(y)] dxdy 



r a/r. 



< Hh^[ ,a]) N / y°- l yK(y) + u(y)} d y 



< Mh-w^Nr*- 1 / y[u n (y) + u(y)] dy. 



As ||u n ||y < L and ||u||y < Q, using (85) we find that 

c(>(x)b{x,y)S(y)[u n (y) -u(y)] dydx < <t>\\ L - [[oM) N{L + Q)r 6 - 1 < 

a/n r 

Now, by changing the order of integration and using Q and H5) we get 

a/n oo 

■«„W.,»>s<»W»>** 



x 
a/n y 



(p(x)b(x, y)S(y)u(y) dx dy 







a/n 



<f>(x)b(x, y)S(y)u(y) dx dy 



j/n 



a/n oo a/n 

< Uh°°([0M) N J A(y)<%+IMU<»([o,a]) J J y x 2a x- 2a b{x,y)u{y)dxdy 

a/n 

a/n oo 
2a / „,0-2a 



< \\<j>h°°([0M) Nn / V Hy)dy+H\\L~([o,a])Cn-^ / y"-^ U {y)dy 



a/n 
2a I „,9-2a 



= n\\ L ™([o,a])( N + C ) n / V°~ 2a <y)dy< ||^|U» ([0 , a] )(JV + C)n-^Q. 



In a similar way we also have 

<p{x)b{x,y)S(y)[u n (y)-u{y)\ dydx 

a/n X 



(86) 



(87) 



< II^IU~([o : a])^ 9 / \u n (y)-u(y)\dy^0 as 
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Now, from (86 )-(88 1 we find that 



< 



cj)(x) [M2(u n )(x) - Mi{u)] dx 
cf>(x) [Ml l (u n )(x) - M 4 (u)} dx 



a/r 



a/n 



<j>(x)M 4 (u)dx 



a/n x 



< 



a/n n 
a r 



(x)b(x , y) S '(y) [u n (x) — u(x)] dy dx + J J <f>(x)b(x , y) S '(y) [u n (x) — u(x)] dy dx 

a / 11 r 

a/n oo 

<>(.r }h(r. ij)S( ij)ii"(.r)dij d.r + J J (j)(x)b(x,y)S(y)u(y) dy dx 

x 



4>(x)b(x , y) S (y) [u n (x) — u(x)} dydx 



a/n x 



+ \\<f>\\L^([0.a]) N ( L + Q) ne 1 + ||0||i-([O,a])(^V + C)?l 2<T Q -> £ RS 11 V OO 

Since <fi and e are arbitrarily chosen, it follows that 

M%{u n ){x) -± M 3 (u)(x) in L^OjoD as n -> oo. 
Case i = 4: By using H4) we have 

\<f>(x)S(x)\ < ||0|U-([o,a])a 8 for a.e.x€[0,a]. 
Then, as <j>(x)S(x) € L°°([0,oo[) for x € [0, a] we find that 



<j)(x) [M2(u n )(x) ~ M 4 (u)(x)} dx 



(89) 



< 



4>(x) [S n (x)u n {x) — S(x)u{x)] dx 



o 

a/n 



[S n {x)u n (x) — S(x)u(x)] dx 



a/n 



4>{x)S{x)u{x) dx 



a /r< 



(j}(x)S(x) [u n {x)-u(x)} dx 



cr/r 



< 



|]0l|i«'([O,a]) CF 71 I u{x)dx + 



(x)S(x) [u n {x)~u(x)} dx 



a/n 



as n — » oo. 



Since (/) is arbitrarily chosen, it follows that 

M2(u n ){x) ^M 4 (u)(x) in L x (]Q,a[) as n -> oo. 



Lemma 3.1 follows from (73), (84l, (891, and (90 1. 



3.2 The existence result 



(90) 
□ 



Theorem 3.2 Suppose that Hl)-H4) hold and assume that uq 6 Y + . TTien f/ie problem |^)-|^|) /ias a solution 
ueCu([0,oo[,Y + ). 

Proof. Choose T, to > 0, and let (u™) ngN be the weakly convergent subsequence of approximating solutions 
obtained in Lemma 2.7 For t e [0, T] we obtain by weak convergence and Lemma |2.6[ i) 

m rn 

/ xu(x,t)dx — lim / xu n (x,t)dx < L(T) < oo, 
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and 



x~ a u(x, t) dx = lim / x~ a u n {x, t) dx < L(T) < oo. 

n— >oo 

1/m 1/r 



Then taking to — > oo implies that it € with ||it||y < 2L(T). Let G L°°(]0, a[). From Lemma 2.7 we have 
for each s <E [0, t] 



u"(s) - 1 u(s) in i 1 (]0,a[) as n — > oo. 



From Lemma 2.7 and Lemma 3.1 for each s £ [0,t] we have 

0) [M"(u"(s))(a;) - M(u(s))(x)] dx ->■ as n -> oo. 
Also, for s G [0, i], using Lemma [2.6[ i) and ||it||y < 2L(T) we find that 
|</»(x)| |M"(u n (s))(a;) - M(«(s))(x)| 



(91) 



(92) 



< H0lU-(]O,a[) 



K(x - y,y) [u n (x - y, s)u n (y,s) + u(x - y, s)u(y, s)] dy dx 



. 

a n—x a oo 

+ J J K (x , y)u n (x , s)u n (y , s) dy dx + J J K(x,y)u(x, s)u(y, s) dy dx 



an a oo 

b( x ,y) s {y) u,l (yi s ) d v dx + J J K x >y) s (y) u (y' s ) dydx 

Ox Ox 



+ / S(x) [u n (x, s) + u(x, s)] dx 



<3| 



» ao ,a D [10L(T) + (N + 1)]L(T). 



(93) 



Since the left-hand side of ( 93 1 is in i 1 (]0,t[) we have by (92), (93) and the dominated convergence theorem 



(j>{x) [M n (u n (s))(x) - M(u(s))(x)] dxds 







— > as n — > oo. 



(94) 



Since <f> was chosen arbitrarily the limit ( 94 1 holds for all 4> £ L°°(]0, of). By Fubini's Theorem we get 
* t 



M n (u n (s))(x)ds ^ M(u(s)){x)ds in L x Q0,a[) as n -> oo. 



(95) 



From the definition of M n for t G [0, T] 
t 

u n (t)= J M n (u n (s))ds + u n (0), 
o 



and thus it follows by (95), (91 1 and the uniqueness of weak limits that 

u{t) = J M(u[s))ds + u(0). 
o 

It follows from the fact that T and a are arbitrary that u is a solution to (|T|) on Cg([0, oof, Y + ). And this 
completes the proof of Theorem |3.2| □ 
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4 Uniqueness of Solutions 

Theorem 4.1 If HI), H2), H3'), H4'), and H5) hold then the problem |7])-|I|) has a unique solution u 6 

c B ([o,oo[,r+). 

Proof: Let us consider u\ and u 2 to be solutions to (JTJ)— (|2j) on [0, T] for T > arbitrarily chosen, with 
Ui(x, 0) = U2(x, 0) and set U = u± — u 2 . We define for n = 1, 2, 3, . . . 



m n (t) = J (x- 17 +x x - a )\U(x,t)\dx. 
a 

Now. for S £ R, we define sgn (<5) as follows 

( 1 <5>0 
sgn (<S) = J 5 = 
{ 1 J < 0. 

Note that if £(•) is an absolutely continuous function of t, then so is 1 1— >■ |C(i)|, an d 

|| C (t)|=sgn(C(t))|cW a.e. (96) 

Also note that by Definition |l.l| iv), u(a;, •) is absolutely continuous on [0, T] for a.e. x € [0, oo[ and therefore 
u(x,t) satisfies Q for a.e. t 6 [0,T], see [9] page 354]. Then, taking the difference of the derivative of the 
solutions u\ and u 2 in we have 



dU(x,t) _ 1 
Ft ~ 2 



- y> 2/) - V, t)ui(y, t) -u 2 (x- y, t)u 2 (y, t)} dy 



K(x,y) [«!(a;,t)«i(j/,t) - u 2 (x, i)ti 2 (y, i)] + / b(x,y)S(y) [u x {y,t) - u 2 (y,t)]dy 



-S(x) [ui(x,t) — 1*2(2;, t)] 



Applying (96 1 we find that 



d\U(x,t)\ 

at 



= sgn(U(x,t)) 



K(x - y,y)[u 1 {x - y,t)u x {y,t) - u 2 (x - y, t)u 2 (y, t)} dy 



K(x,y) [u\{x,t)ux{y,t) - u 2 (x,i)u 2 (y,t)]dy 



+ / b(x,y)S(y) [u x {y,t) - u 2 (y,t)]dy - S(x) [u^x.t) - u 2 (x,t)] 



Multiplying by (a; CT + x x CT ) and integrating from to t and from to n w.r.t. r and x respectively, and 
applying Fubini's Theorem, for each n and < t < T we obtain 



t n 

m n {t)= J J {x- a + x x -°) sgn (U(x, r)) 




K{x - y,y) [mix ~ t/,r)ui(y,r) - u 2 (x - y,T)u 2 (y,T)} dy 



K{x,y) [ui(x,T)ui(y,r) - u 2 (x,T)u 2 (y,r)] dy 



+ / K x ^y) s {v) [ui(y,t) - u 2 {y,t)] dy - S(x) [ux{x,t) - u 2 (x,t)] 



dx dr. 



(97) 
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Using the substitution y — x = x' in the first of the inner integrals w.r.t. x and y on the right hand side of (97 1 
we find that it becomes 

n x 

{x~ a + x A ~ CT )sgn (U(x,t))^ J K(x - y,y) [u\{x - y,T)u\{y,T) - u 2 (x - y, t)u 2 (y, r)] dydx 
o o 

n n—x 

- [(x + yy a + (x + y) x - a ] sga(U(x + y,r))K(x, y) [ux(x,r)ui(y,T) - u 2 (x,T)u 2 (y,T)]dydx. 

o o 



Inserting this into ( 97 1 gives 

[(x + y)-° + (x + y) A - ff ] sgn (U(x + y, r)) - {x~° + x^) sgn (*7(z, r)) 



t n n—x 

'l 



m"(i) = 



■K(x,y) [ut(x,T)ui(y,r) - u 2 (x,T)u 2 (y,r)}dy dx dr 

t n oo 

(x~ a + x A_cr ) sgn (t/(#, r))K{x,y) [ui(x, r)ui(y, r) - w 2 (^, t)u 2 (]J, t)] dydxdr 

n-a; 
£ n oo 

(x _,T + x A ~ CT ) sgn (U(x,T))b(x,y)S(y) [tti(y,r) - «2(2/,t)] dydxdr 

a o x 

t n 

( x - a +x x - a )sgn(U(x,T))S(x) [ui(x,r)-u 2 (x,T)]dxdT. (98) 



We note now. using the symmetry of X. by interchanging the order of integration and interchanging the roles 
of x and y yields the identity 

n n—x 

(x~ a + x x ' a ) sgn (U(x,r))K(x,y) [ui(x, r)ui(y, r) - u 2 (x,r)u 2 (y,r)] dydx 

o o 

n n—x 

{if" + V X ~ a ) sgn (U(y, t))K(x, y) [ui(x, r)ui(y, r) - u 2 (x, r)u 2 (y, r)] dydx. (99) 
o o 

For a;, y > and i S [0, T] we define the function w by 

™(x, y, t) = [{x + y)- a + (x + y) x -°] sgn (U (x + y, t)) 

-{x-° + x x -°) sgn (U(x, t)) - (y- a + y x -°) sgn (U(y, t)) . 



Using (99 1 and this definition, we can rewrite (98 1 as 

t n n—x 

TO ™W=o / / / w ( x 'V' T ) K ( x >y) i u i( x , T W(y,T) + u 2 (y,T)U(x,T)] dydxdr 



2 





t n oo 

{x~ a + x x - a ) sgn (U(x, r))b(x, y)S(y)U(y, r)dy dx dr 

x 
t n 

{x~ a + x x ~°) sgn (U(x, t))S{x)U(x, T )dxd T 



t n oo 

X-c 



(x- a + x x - a ) sgn (U(x, t))K(x, y) [ui(s, r)U(y, r) + u 2 (y, t)U(x, t)] dy dx dr. (100) 



n—x 
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Since the third integral and the second term in the fourth integral of (100 1 are positive, we can delete it and 
get the following estimate 



t n n—x 



m n (t) < 1 
w - 2 





t n n—x 



+ - 



+ 





t n oo 



x 
t n oo 



w(x, y, t)K{x, y)u\(x, r)U(y, r) dy dx dr 

w(x, y, t)K(x, y)u2(y, t)U(x, t) dy dx dr 
gn (U(x, r))b(x, y)S(y)U(y, r)dy dx dr 

+ x x ~ a ) sgn (U (x, t))K(x, y)u\(x, r)U(y, r) dy dx dr 



o n-x 



= / [hi(r) + h%{r) + 7 33 (r) + J 34 (r)] dr. 



(101) 



Taking into account that for all q € K \q\ — qsgn(q) holds, we can estimate 

w{x, y, t)U(y, t) <[[(x + y)-° + (x + y) x -] + (x~" + x x ^) - (y- + y x -)] \U(y, t)\ 



Using the inequalities (|9| and (11 1 we find that 

w(x,y,t)U(y,t) < [[2— \ x - + y-) + ( x x - + y x -)] + { x - +3 *-°) _ ( y - + y x -)] \U(y,t)\ 

< [[(x-° + y-°) + (x x -°+y x -°)]+(x-°+x x -°)-(y- a + V X - a )] \U(y,t)\ 

< 2(x-° + x x -°)\U(y,t)\. (102) 



Now, we use ( 102 1 to work on each term of the right hand side of ( 101 1 



t n n—x 



hi(T)dT< / / / (x-° + x x - a )K{x 1 y)u 1 (x,T)\U(y,T)\dydxdT. 



(103) 



o 



Using the estimate H3') for K(x,y) and (10 1 we get 

t t n n—x 



Jhi(r)dT< Kl J J J {x~° +x x -°){x-° +x x -< T )(y-' T + y X - a )u 1 {x,T)\U(y,r)\dydxdT 

{x-° + x x -°) 2 u 1 (x,T)(y- a + y x -°)\U{y,T)\dydxdT 





t n n—x 





t n n—x 

<2 Kl J J J (V 2CT + x 2{ - x -^) Ul {x,T)(y-° +y x -°)\U(y,T)\dydxdT. 



Due to A — a, a G [0, 1/2] and the definition of m n (t) from the inequality above it follows that 
t t 



h\ij)dT < 2ki / m n (r) 



(x- 2a +X 2{X - a ^U 1 {x,T)dx+ I (x- 2a +X 2{X - a) ^jU 1 {x i T)dx 



dr 



< 2ki / m n (r) 



(x 1 + l)ui(x, t) dx + I (1 + x)ui(x, t) dx 



dr 



< 2ki / m n (r) 



2 x 1 ui(x,r)dx + 2 / xu\{x, t) dx 



dr<Ai J m"(r) dr, 



(104) 



where Ai = Ak\ sup se [ t i ||ui(s)||y. In the same way, there is a constant A 2 such that 
t t 
I 32 (t) dr < A 2 J m"(r)dr. 
o o 
Now, changing the order of integration in J33 we have 

t t n 00 

h?,{T)dT= I I sga.(U(x,T))b(x,y)S(y)U(y,T)dydxdr 



X 



< 



< 



n y 00 n 

K x > y) s {y)P(y, r)\dx dy + 

n 

00 y 

K x ^y) s {y)\u{y,T)\dxdy- 

n 





n J/ 



6 ( a; ^)'S'(y)|C / (y I ' r )|rfa;c ; y 

K x ^y) s {y)\u{y^ T )\dxdy 



dr 



dr. 



By using p) and H4') we find that 



L 



/ h 3 (r)dT < J N J y x -°\U{y,T)\dy + N J yy^piy^dx 

n 

00 

Nm n { T ) + Nn 9 - 1 J y[u x {y, r) + u 2 (y, r)]dy 
< I [Nm n {T) + Nri 6 ' 1 + \\u 2 \\ Y )] dr 



dr 



< 



d,T 



N J m n {r)dT + Nrf-' (||ui||y + ||u 2 ||y)i. 




-34 
00 00 







To consider J34 we first see that 

(aT CT + x x - a ) sgn (U(x, t))K(x, y)ui{x, t)U(y, t) dy dx 

OO OO 

<2 Kl J J (x~ 2 ° + x 2 ^) (y-° + y x -°)u 1 (x,t)\U(y,t)\dydx < 


Thus, the dominated convergence theorem leads to 
J hi( T ) dr — !• as n — >• 00. 



0G. 



Therefore, due to ( 101 1, ( 104 !-( 107 1 and taking A = Ax + A 2 + N we obtain 



m(t) := / (x- a + x x ~ a ) \U(x,t)\dx = lim m n (t) 

J n— >oo 




< lim / [7 31 (r) + J 32 (r) + 7 33 (t) + 7 34 (t)] dr 



t t 

< lim A / m™(r)dr+ lim / 7 34 (r) dr 



f 00 




(105) 



(106) 



(107) 
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From where we have the inequality 
t 

m(t) < A J m(r) dr. (108) 
o 



Applying Gronwall's inequality, see e.g. (l4[ page 361], we obtain 

oo 

m(t) = I (x- a + x x - a ) \U(x,t)\dx = for all t € [0,T]. 



o 

Thus, we have that 

ui(x,t) = U2(x,t) for a.e. x€]0,oo[. 

□ 
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